Splitting for Dissipative Particle Dynamics

We study numerical methods for dissipative particle dynamics, a system of stochastic differential equations for simulating particles interacting pairwise according to a soft potential at constant temperature where the total momentum is conserved. We introduce splitting methods and examine the behavior of these methods experimentally. The performance of the methods, particularly temperature control, is compared to the modified velocity Verlet method used in many previous papers.